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Extensive Air Showers are still the only source of information on primary cosmic rays and their interactions 
at energies above PeV. However, this information is hidden inside the multiplicative character of the cascading 
process. Inspite of the great experimental and theoretical efforts the results of different studies are often ambiguous 
and even conflicting. These controversies can partly be referred to imperfections of our models of high energy 
particle interactions. 

The first part of the paper is concerned with this problem. The author thinks that the present models should be 
corrected to give slightly deeper penetration of the cascade into the atmosphere. In this respect the modification 
suggested by the QGSJET- II model seems to be the step in the right direction. The Sibyll 2.1 model provides 
a similar penetrating properies. However, this modification is not enough and a small additional transfer of the 
energy from EAS hadrons to the electromagnetic component is needed too. As a possible candidate for such a 
process the inelastic charge exchange of pions is discussed. 

In the second part of the paper the author discusses the need to account for the interaction of EAS with the 
stuff of detectors, their environment and the ground in the light of the 'neutron thunder' phenomenon, discovered 
recently. 



1. Introduction 

There is a big progress in the analysis of ex- 
perimental data on extensive air showers ( EAS ) 
during the last two decades. However, one cannot 
say that we understand all the phenomena and 
characteristics of EAS which we observe. Partly 
this dissatisfaction is due to the controversies in 
experimental data themselves, partly due to still 
remaining imperfections of the analysis. We cer- 
tainly need to improve our understanding of EAS. 

This paper does not aim to give a comprehen- 
sive review of all high energy interaction models, 
event generators and EAS simulation codes. It 
consists of two different parts. In the first part I 
point out some problems related to the particle 
interaction models which so far pose questions 
at high energies. I do not go into the theoreti- 
cal foundations of various interaction models, but 
stay within a pure phenomcnological approach. 
Within it I indicate the possible way to improve 
the models. The theoretical basis of some recent 
models can be found in pQ. 

In the second part of the paper I shall touch 
the problems related to some effects of the EAS 
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interaction with the environment. 

2. High energy interactions 

2.1. The consistency of the results 

The EAS is a complex phenomenon - it has 
several different components: electromagnetic - 
electrons, positrons and gamma-quanta, muons 
and hadrons - nucleons, pions, kaons and so on. 
Besides that there are neutrinos which need mas- 
sive detectors to be studied. Due to their small 
interaction cross-section they are detected not as 
multiple shower neutrinos, but as single ones. So 
far they are not combined with other EAS com- 
ponents in the analysis of experimental data, but 
they certainly play a role in the energy balance. 
Optical cherenkov and fluorescence photons emit- 
ted by charged shower particles are also used as 
a powerful tool for the study. 

Since the characteristics of observed showers 
are the product of the primary cosmic-ray ( CR ) 
energy spectrum, mass composition and high 
energy interactions, the only way to disentan- 
gle them is to achieve the self-consistency in 
the derivation of the properties of primary CR 
from different observables and vice versa - the 
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derivation of observed characteristics for differ- 
ent shower components and different observation 
levels from the same primary CR and the inter- 
action model. 

There were many efforts in the past to use mod- 
els of the popular CKP- or scaling type. With the 
development of the QGS model [2] it has been 
shown that this model gives a satisfactory de- 
scription of both EAS [314] and single, unasso- 
ciated CR components in the atmosphere [51617] , 
However, those old studies used as a rule different 
cascading algorithms and programs, which cer- 
tainly produced an additional uncertainty in the 
results and reduced the credibility of the conclu- 
sions. It is to the credit of the KASCADE people 
who spend great efforts to develop and to dis- 
tribute freely the CORSIKA code [8]. With this 
code the analysis of experimental data can now 
be made at the level much better than before. 

2.2. The improvement of models 

An early analysis of models indicated that the 
best consistency for the mean logarithmic mass 
(InA) of primary CR derived from the jf- ratio 
and from X max can be achieved for the QGSJET 
model [5]. Here N^,N e are muon and electron 
sizes of EAS respectively and X max - the depth 
of the shower maximum. Later a similar conclu- 
sion about the preference of QGSJET model has 
been made on the basis of the analysis of the EAS 
hadronic core [10] . After some improvements the 
SIBYLL model, version 2.1 joined the list of the 
best, most popular and often used models [TTj . 

However, the closer look reveals that some in- 
consistencies still remain. It has to be said that 
indications of possible inconsistencies appeared 
more that 30 years ago when the mismatch be- 
tween the direct and indirect measurements of the 
primary energy spectrum has been noticed: the 
indirect measurements based on the EAS model 
calculations gave as a rule the higher CR inten- 
sity in the PeV region than that derived by the 
extrapolation of direct measurements from the 
lower energies - the so called 'bump' problem [12] . 
More recently this mismatch has been confirmed 
by [13]. Among possible explanations there was 
an assumption that even the best models give an 
overestimation of the primary energy from the ob- 



servations in the atmosphere. It could happen if 
the shower penetrates deeper into the atmosphere 
and has more charged particles at the observation 
level than it is expected from model calculations. 

Observations of the EAS cherenkov light in 
the PeV region confirmed this deeper penetra- 
tion [14] . As a consequence, the primary mass 
attributed to such showers derived from observed 
X max values and N^/Ne ratio after the compari- 
son with model calculations turned to be smaller 
than the true primary mass. There was a num- 
ber of ideas how to increase the penetrability of 
the showers, for instance, introducing the higher 
cross-section for the charm production [15] or 
hypothetical strangelets [16117] . but those mod- 
els are still in the stage of development. The 
possibility to improve the models were discussed 
also in [14|18] , In [14] it has been assumed that 
the cross-section and the inelasticity of the pro- 
ton interactions in the air are in fact smaller 
than in the models, although they still agree 
with measurements at the lower end of the er- 
ror bars. Their reduction allowed to improve the 
agreement between the predictions of the mod- 
els and the results of the X max measurements. 
There were some indications of the lower cross- 
sections in the past measurements of hadrons in 
the EAS cores [19]. The latest measurements of 
the inelastic cross-sections confirmed the slower 
rise of the interaction cross-section with energy 
20 21 . Therefore, there are experimental indi- 
cations that EAS may in fact penetrate deeper, 
than predicted by models. 

There are also efforts to improve models not 
just on the pure phenomenological, but also on 
the theoretical basis. The idea that the density 
of partons at high energies becomes so high that 
they cannot interact independently of each other 
has been discussed long ago [22]. However, it is 
to the credit of S.S.Ostapchenko, who updated 
the QGSJET01 model including the non-linear 
effects of parton interactions, developed it to the 
status of the Monte Carlo event generator and to- 
gether with his colleagues in Karlsruhe incorpo- 
rated it into the Corsika code [23|24] , As a con- 
sequence of the non-linear effects, the interaction 
cross-section ( at least for pions ) , the multiplicity 
of secondaries and the inelasticity of the collisions 
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decreased slightly which helped atmospheric cas- 
cades to penetrate deeper. Apparently the reduc- 
tion of the inelasticity plays the major role in the 
increased penetrability. Due to its smaller inelas- 
ticity the updated Sibyll 2.1 model also provides 
EAS with a greater penetrability than previous 
models. Certainly these improvements are the 
step in the right direction. 

However, the only introduction of the non- 
linear effects of parton interactions seems to be 
not enough. This suspicion appears when the ex- 
amination of the hadron component is included 
into the analysis. It has been shown in [25] that 
the primary mass composition derived mainly 
from hadron and muon components is heavier 
than that which can be obtained using mainly 
electromagnetic and muon components. Muons 
are usually less model dependent at the fixed pri- 
mary energy, since they are penetrating particles 
and are collected from all atmospheric altitudes 
representing something like an integral over the 
longitudinal profile of the shower. Taking them 
as the basis for the comparison we should expect 
that for well tuned, consistent models the analysis 
of ratios Ne/N^ and Nh/N^ should give the same 
(InA). The larger (InA) value (2.25±0.08) found 
in KASCADE experiment from Nh/N^ analysis 
than that from N e /N^ (1.90±0.05) [25] indicates 
that the difference between N e /N^ and Nh/N^ 
in the present models is too low ( Figure 1 ). In 
the more realistic model this difference should be 
increased. 

The theoretical basis of QGSJET-II model be- 
sides the reduction of the interaction cross sec- 
tion for pions and an inelasticity requires also 
the reduction of the multiplicity of the secondary 
particles. Simulations of EAS showed that tran- 
sition from QGSJETOl to QGSJET-II model at 
the fixed energy leads to the rise of Ne/N^ and 
Nh/N^ ratios both for primary protons and for all 
primary nuclei. It is because the EAS electromag- 
netic and hadron components follow each other 
in the lower part of the atmosphere, i.e. beyond 
the shower maximum, in an approximate equilib- 
rium. However, the difference between {InA) val- 
ues derived from N e /N^ and N^/N^ ratios in the 
QGSJET-II model becomes even larger than for 
QGSJETOl. While the (InA) value derived from 
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Fi gure 1 . The illustration of the way to get the consistent 
value of mean logarithmic mass (InA) from (a) -/V e /7V M and 
(b) Nh/N^ measurements. Initial values of TVe/iV^ and 
Nh/N^ are calculated for 1 PeV primary protons, nitrogen 
and iron nuclei using the QGSJETOl model and shown by 
the full line. The dashed line shows the same ratios for 
the QGSJET-II model: N^/N^ and N h /N^ rise by the 
same factor. The dash-dotted line demonstates the effect 
of the charge exchange combined with the QGSJET-II: 
Ne/Nfj, rises, but N^/N^ falls. The thick arrow indicates 
the consistent value of (InA). 



Ne/Np ratio rises from 1.9 to 2.52, that derived 
from Nh/Nfn rises from 2.25 to 3.47 ( see Figure 
1 ). Therefore, the reduction of the interaction 
cross-section for pions, of the inelasticity and the 
multiplicity in the QGSJET-II, does not remove 
the existing difference between (InA) values and 
the inconsistency still holds. 

To eliminate this inconsistency we have sug- 
gested the additional transfer of the energy into 
an electromagnetic component in the cascading 
process [18]. It has been made on a pure phe- 
nomenological basis. Later we have suggested the 
so called 'sling effect' in nucleus-nucleus interac- 
tions as the process responsible for a deeper pen- 
etration into the atmosphere of cascades induced 
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by primary nuclei and also for an additional elec- 
tromagnetic component in them [26127] . This ef- 
fect could serve as a possible theoretical basis of 
the needed model modifications. However, this 
effect being very probable, seems to be small to 
give the noticable changes. The needed effect has 
to be stronger and we now think that the charge 
exchange process is the likely culprit. 

In fact charge exchange processes are already 
taken into account in both QGSJET01 and 
QGSJET-II models. The question is whether it is 
possible to modify the probability of this process 
without a conflict with the existing experimental 
data and whether such a modification could give 
the consistent (InA) value both from N e /N^ and 
Nh/Np ratios. 

In order to analyse this possibility I used 
the option provided by the HDPM model. It 
is the only model within the CORSIKA code 
where one can switch on and off the charge ex- 
change processes and by this way to estimate 
the effect which this process has on N e /N^ 
and Nh/N^ ratios. Actually I determined the 
ratios R elx = {N c /N^)hdpm + / (N e / N ^)hdpm_ 
and R htl = {N h /N t _ l ) HD p M+ /{Nh/N^hdpm- 
for the HDPM model with and without 
the charge exchange process ( denoted as 
HDPM+ and HDPM_ respectively ). Then 
I applied these ratios to the QGSJET-II 
model as log^^/N^QGSJET-n+ch.exch. = 
log^Ne^/N^QGsjET-n + elogRe.h- Here the co- 
efficient e limited by < e < 1 gives an estimate 
of possible increase of the charge exchange needed 
to get a consistent value of (InA) derived from 
N e /N^ and Nt/N^ ratios, denoted as (lnA) e)1 
and (InA) hp, respectively. 

This exersise shows that by this way it is pos- 
sible to achieve the needed consistency. Sim- 
ulations show that inspite of the addtionsl en- 
ergy transfer the shower size at the sea level re- 
mains practically the same as in the absence of 
the charge exchange process, whereas the num- 
ber of hadrons and muons decreased. So that 
in contrast to the modifications provided only by 
non-linear effects in the QGSJET-II model when 
both Ne/Nft and Nh/N^ rise, the increase of the 
charge exchange probability leads to the rise of 
Ne/Np, but to the fall of Nh/N^, because both 



and Nh fall but the latter falls stronger. For 
example, with e — 1 (lnA) e ^ rises from 2.52 to 
2.98, but (lnA)hn falls from 3.47 down to 2.86 
( see Figure 1 ). This 'overshooting' is due to that 
the charge exchange process is already taken into 
account in the QGSJET-II model and applica- 
tion of the expression given in the previous para- 
graph with e = 1 makes this process too strong. 
The consistent value of (lnA) eil — (lnA) efl = 
2.94 ± 0.09 is achieved at e = 0.88 ± 0.12. The 
errors are derived from the statistical errors of 
mean values of iV e /iV M and Nh/N^ ratios obtained 
by Monte-Carlo simulations with CORSIKA6.014 
and CORSIKA6500 codes and taking the values 
of (lnA) eii = 1.90±0.05 and (lnA) hf , = 2.25±0.08 
obtained from the experimental data using the 
QGSJET01 model 25J. Systematic errors are dif- 
ficult to evaluate at this stage of analysis, but ac- 
tually our estimates are given just to demonstrate 
the principal opportunity to use the QGSJET-II 
model with an enhanced charge exchange proba- 
bility to get a consistent estimates of the primary 
CR mass composition. 

The higher consistent value of (InA) = 
2.94±0.09 compared with values of 1.90±0.05 and 
2.25±0.08 derived from A e /A M and N h /N^ ratios 
with the QGSJET01 model means that the true 
primary CR mass composition should be heavier 
in the new analysis. 

Simulations show also that with the increased 
charge exchange the depth of maximum for pro- 
ton induced showers shifts upwards by about 5 
gcm~ 2 , but for nitrogen and iron induced show- 
ers it moves downwards by 1 gem -2 and 5 gem -2 
respectively. Hence the increase of the charge 
exchange cannot destroy the deeper penetration 
of EAS by about 20 gem" 2 provided by the 
QGSJET-II model. 

2.3. The inelastic charge exchange for pi- 
ons 

Now we shall discuss the process which could 
be used to increase the charge exchange probabil- 
ity. At the end of sixties one of the inventors of 
the ionization calorimeter, V.S.Murzin from the 
Moscow University used this detector to study in- 
teractions of CR pions. He has found that with 
a considerable probability the charged pion could 
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lose its charge but preserve a good part of its ini- 
tial energy. He called this process 'the inelastic 
charge exchange' [28] . Actual numbers were the 
following: the probability to preserve more than 
0.5 of the energy in the collision of the pion with 
iron nuclei has been estimated as 10% and this 
probability seemed to be independent of energy 
in hundred GeV - TeV energy range. 

Since then numerous experiments have been 
made on the production of neutral pions in 
hadron interactions. Their results can be 
found in the Particle Data Group archive 
( |http://durpdg.dur.ac.uk/HEPDATA/reac.html| 
( see also the list of literature in [29] ). However, 
measurements of collisions with nuclei are still 
sparse and cover mostly the high transverse mo- 
mentum ( P t ) region. The relevant data relate 
mostly to n^P interactions and often give the val- 
ues of invariant cross-sections |29|30|31|32|33|34] . 
After integrating over P t some extreme examples 
30 31 33J of the obtained inclusive spectra for 
neutral pions are shown in Figure 2. 

It is seen that the results of different experi- 
ments have a considerable spread. While the re- 
sult of one experiment matches the model per- 
fectly [33j . the 7r°-spectrum obtained in another 
experiment follows the spectrum of n + [30] and 
there is a spectrum which is definitely higher than 
that which matches the model [31] although they 
agree with each other within error bars. This 
difference cannot be caused by the energy de- 
pendence of the process since models and experi- 
ments show no appreciable dependence in the 58 - 
360 GeV energy interval [34] . The examination of 
these data indicates that the probability to have 
7T° with X > 0.5 is not 10%, but about 7% at 
these energies and for X > 0.7 it is about 2.3%. 

It seems that on the basis of the spread of ex- 
perimental data and the absence of the preference 
between the different results the model can be re- 
tuned according to indications of the EAS anal- 
ysis towards the higher probability of the charge 
exchange process for pions. This can improve the 
consistency between (InA) values obtained from 
different EAS components. Since at X > 0.7 the 
difference between different experiments rises up 
to ~ 2.5 and continues growing with X, the big 
value for the estimate of e = 0.88 ± 0.12, pre- 
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Figure 2. The inclusive spectrum of neutral pions pro- 
duced in 7r~P interactions. Full circles - for the 58 GeV 
7r~ beam energy [311 . open and full stars - for 80 and 140 
GeV respectively |33l . triangles - for 360 GeV [30]. Hor- 
izontal error bars indicate the interval of X used for the 
determination of the invarint cross-section. Histograms 
are spectra of 7r + ,7r° and n~ calculated with CORSIKA 
INTTEST version for the QGSJET01+GHEISHA interac- 
tion model at 80 GeV. The use of the QGSJET01 model at 
this energy is justified since non-linear effects introduced 
in QGSJET-II appear at much higher energies. 



sented in the previous subsection, is not inconsis- 
tent with these experimental data. 

There is another point which should be re- 
marked. The appearance of leading neutral pi- 
ons in charged pion collisions observed in cosmic- 
ray and in some of accelerator experiments is 
interpreted within the framework of the triple- 
region description with a substantial contribution 
of RRP-term. Its contribution should give the flat 
behaviour of the inclusive 7r° cross-section at large 
X > 0.7. It is not seen in the data shown in Fig- 
ure 2. This discrepancy is not clear. However, in 
can be that 7r°'s with very large X are biased 
in accelerator experiments by triggering condi- 
tions, with partially suppressed low multiplicity 
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and diffraction events as has been mentioned in 
[30] . while the introduced corrections are model 
dependent. There was no such bias in cosmic-ray 
experiments. 

The experimental value of the mean fraction a 7 
of energy transferred by 7r~ to ir° is independent 
of energy and equal to 0.25±0.01 28J. The a 7 
provided by the QGSJET01 model shows a slight 
decrease with the energy and above 10 TeV falls 
below the experimental value ( Figure 3 ). Inde- 
pendently of the origin of this fall ( e.g. change 
of the mass composition of secondaries ) it is an- 
other indication that the energy transfer into the 
electromagnetic component should be slightly in- 
creased to improve the model. 
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Figure 3. The mean fraction of energy a 7 transferred 
to neutral pions in 7r~ P collisions as a function of energy. 
The experimental value of a 7 = 0.25±0.01 is from |28| and 
denoted by full and dashed lines. The full circles show a 7 
values calculated with the QGSJET01+GHEISHA model. 
The calculated values are slightly smaller than experimen- 
tal ones above 10 TeV. 



3. Interactions of EAS with detectors and 
an environment 

The second part of this paper relates to inter- 
actions of EAS particles not with air nuclei, but 
with another target: the stuff of the detectors, 
their environment and the ground. The aim of 
this part is to draw an attention to the possi- 
ble contribution of low energy neutrons created 
in such interactions to the signal in the hydro- 
gen containing detectors, deployed particularly at 
mountains. 

It is well known that the thickness of the shower 
disk depends on the distance from the shower core 
rising from 2-3 m at the EAS center to tens of m 
at about 1 km from it. It corresponds to the 
time 'thickness' of a few hundreds ns. The de- 
layed particles which appear a few microseconds 
after the main shower front were observed and 
discussed long ago [35|36I37j . However, the dis- 
covery of neutrons delayed by hundreds microsec- 
onds in the shower core made by Chubenko A. P. 
and his group with the Tien-Shan neutron moni- 
tor [38|39|40l41|42j seemed to be unexpected and 
attracted an attention |43|44|45|46"|47] . 

Though the dispute on the interpretation of 
this finding is not finished the majority of par- 
ticipants is inclined to explain it by the interac- 
tion of hadrons in the EAS core with the stuff 
of the detector, i.e. by multiplicative processes 
in the lead of the neutron monitor with the sub- 
sequent long diffusion and thermalization of re- 
leased neutrons in the monitor's moderator and 
reflector |47I48|49] ( see, however, another view 
in [50] ). Independently of the interpretation the 
phenomenon is very spectacular and looks as the 
neutron 'thunder' which appears with a time de- 
lay after the 'lightning' which is the EAS itself 
pIBTj . 

If the interpretation of the majority is correct, 
the observed neutrons can be produced also in 
the ground which is not so heavy as lead, but 
nevertheless there are many heavy elements in it 
( mainly Si ) which could produce neutrons be- 
ing hit by an EAS core. There is also water in 
it which serves as a good moderator like in nu- 
clear reactors. The influence of the ground and 
ground-based environment has to be more sub- 
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stantial at the mountain level, where EAS cores 
are much more energetic and a good part of the 
year the ground is covered by snow. Sometimes 
this snow is of meters thick ( Tien-Shan, Aragats, 
Chacaltaya, South Pole etc. ). As for the Tien- 
Shan station an additional factor is that it is built 
on the permafrost with a good part of ice in it. 
Since neutrons can diffuse up to long distances 
from the place where they are produced [52] their 
effect might be noticable even at shallow depths 
underground. 

Many running EAS arrays use water or ice 
cherenkov detectors: Pierre Auger Observatory, 
Milagro, Nevod, Ice- Top etc. At the first 
sight they should not be sensitive to neutrons, 
since they are neutral particles and mostly non- 
relativistic. However, the experimental study 
of the neutron 'thunder' revealed that delayed 
neutrons are accompanied by gamma-quanta and 
electrons |41| , which in principle could give a sig- 
nal in water tanks. 

It is particularly relevant to Pierre Auger Ob- 
servatory. The comprehensive modelling of the 
effect of albedo neutrons emitted by the ground 
as a result of the EAS interaction is complicated 
and planned for the future paper. Here I show 
that the effect can be noticable even taking into 
account only EAS neutrons. 

Simulations of the EeV proton induced show- 
ers observed at the altitude of 1400 m a.s.l. show 
that neutrons are the most abundant among EAS 
hadrons and their lateral distribution function 
( LDF ) is wider than LDF for protons and pi- 
ons. At the typical distance of ~1 km from the 
core the density of neutrons with energy above 50 
MeV and their energy density is about the same 
as for muons of this energy and about 10% of the 
gamma-quanta plus electrons with energy above 
1 MeV ( Figure 4 ). 

Moreover after about 5 [is behind the EAS 
front the neutron component at 1 km from the 
EAS axis becomes dominant, overtaking muons, 
electrons and gamma-quanta ( Figure 5 ). 

This distance and the time delay are right the 
working distances and times for Pierre Auger Ob- 
servatory, so that the possible contribution of 
neutrons to the signal from their water tanks 
should be analysed and taken into account if nec- 
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Figure 4. Lateral distribution of particle numbers (a) 
and the particles energy (b) for 1 EeV primary proton in- 
cident vertically at the level of 1400 m a.s.l. It is seen 
that neutrons could contribute up to 10% to the signal 
of water tanks at 1km distance from the shower axis, if 
among products of their interaction with water are rela- 
tivistic electrons. 



essary. The same remarks could be referred to 
hydrogen containing plastic scintillators used in 
many other large arrays ( Yakutsk, Telescope Ar- 
ray etc. ). As it has been mentioned above, signals 
delayed by /is ( 'subluminal pulses' ) have been 
already observed in large scintillator arrays, such 
as Volcano Ranch [36|37j . 

A good analysis of the possible effect of de- 
layed particles on the primary energy estimation 
has been made in 53J applicable to the AGASA 
array. It has been shown that overestimate of the 
primary energy for its scintillators and the acqui- 
sition system cannot exceed a few percents. For 
other arrays it may be higher. 

Presumably the effect of 'the neutron thunder' 
can be applied in practice for the neutron caro- 
tage of the upper layers of the ground. Instead 
of the artificial neutron source in this method the 
ordinary EAS can be used since EAS cores carry 
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Figure 5. Arrival time distribution of electromagnetic, 
muon and neutron component of the shower at core dis- 
tances less than 10m (a), 100m (b) and 1000m (c). It 
is seen that at 1000m from the core neutrons dominate 
among other particles after 5 p.s. 



on and produce a lot of secondary neutrons. Also 
'the neutron thunder' can be used for the search 
of water on the Moon [54] or on the surface of 
other planets |55|56j . 

4. Conclusion 

The analysis of existing controversies in the in- 
terpretation of experimental data on EAS indi- 
cates that an improvement of our understanding 
of the EAS phenomenon and the self-consistency 
of results on primary CR derived from EAS can 
be achieved by a moderate modification of the 
current particle interaction models. This mod- 
ification has to result in a slightly deeper pen- 
etration of EAS into the atmosphere as well as 
in the increased transfer of the energy from the 
hadronic to electromagnetic components of EAS. 
The account for non-linear effects in parton in- 
teractions like that in the QGSJET-II model and 
an increasing probability of inelastic charge ex- 



change processes for pions can help. 

Here it is appropriate to make some general re- 
marks. The nuclear and electromagnetic nature 
of EAS has been established at the end of forties. 
That was the time when the world greatest accel- 
erator - Dubna Synchrophasotron had not been 
commissioned and CR were the unique source of 
information about high energy interactions. It is 
surprising that now, after about 50 years of the 
leading role of accelerators in the field, CR are 
still able to contribute to our understanding of 
the high energy interactions. Another point is 
that after nearly 70 years since the discovery of 
EAS by Pierre Auger and Roland Maze we still 
develop our understanding of this phenomenon. 
The discovery of the 'neutron thunder' certainly 
complements our knowledge of the EAS develop- 
ment and is worth of further experimental and 
theoretical study. 
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